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Abstract: Adomian polynomials (AP's) are expressed in terms of new objects 
called reduced polynomials (RP's). These new objects, which carry two subscripts, 
are independent of the form of the nonlinear operator. Apart from the well-known 
two properties of AP's, curiously enough no further properties are discussed in the 
literature. We derive and discuss in full detail the properties of the RP's and AP's. 
We focus on the case where the nonlinear operator depends on one variable and 
construct the most general analytical expressions of the RP's for small values of 
the difference of their subscripts. It is shown that each RP depends on a number 
of functions equal to the difference of its subscripts plus one. These new proper- 
ties lead to implement a dramatically simple and compact Mathematica program 
for the derivation of individual RP's and AP's in their general forms and provide 
useful hints for elegant hand calculations of AP's. Application of the program is 
considered. 

Keywords: Adomian decomposition method, Adomian polynomials, diophan- 
tine equations, Mathematica. 



1 Introduction 



Decomposition methods continue to develop and gain ground in applied mathemat- 
ics and integral methods. They have been reviewed, modified and applied to differ- 
ent fields of science and engineering. The idea behind a decomposition method, 
which is used to solve differential, integro-differential, algebraic equations and 
their combinations, is to obtain by easily handled successive iterations or recur- 
sions approximate solutions within a predefined accuracy. From this point of view, 
the most applied methods are a variety of domain decomposition methods, which 
solve boundary value problems [Davies, Crann, Kane, and Lai~| ( |2007| ); |Han, and 



Atluri <|2007l; Huang, Zhang, Ma, and Wang (2008) 


Patricio, Mattheij, and de 


With (2008 1; Takei, Yoshimura, and Kanayama (2008 1; 


Vodicka, Mantic, and Paris 


1 Baskent University, Department of Mathematics, Baglica Campus, Ankara, Turkey. 
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[ (2007 1], and Adomian decomposition method and its modifications, which tackle 



both initial and boundary values problems | 


Adomian ( 1994); Adomian, and Rach 


(1992) 


; Cochran ( 1968 1; El-Tawil, Bahnasawi, and Abdel-Naby |(|2004|);|Hosseini 


(2006 1; Lai, Chen, and Hsu 


(2008); 


Wazwaz 


(2000)); Wazwaz and El-Sayed 


(2001) 


; Zhang (2005)]. 



Adomian decomposition method and some of its modifications and extensions are 
mathematical tools providing analytical and rapidly convergent (if not exact) solu- 
tions to a variety of problems in nonlinear science. A problem in nonlinear science 
is modeled by 



Lu + Ru + Nu = g , 



(1) 



where the unknown (scalar, vector or matrix) function u(x) is subject to some initial 
and/or boundary conditions. The function g(x) is a source term, (L + R) is a linear 
operator and N includes nonlinear operator terms, if there are any. The original 
idea due to Adomian is to split the linear operator (L + 7?) into two linear terms, L 
and R, where L represents the highest order derivative and is easily invertible and R 
groups the remaining lower order derivatives. 

Integrating ([T]) leads to a Volterra integral equation [Cochran (1968)] 



f-L- l (Ru)-L- l (Nu 



(2) 



where / is the sum of L 1 (g) and the terms arising from the application of the initial 
and/or boundary conditions to u. It seems that a part of the method originated in 



the work of Cochran | Cochran ( 1968 1]. The other new part of the method consists 
in replacing u by sum of components u m 



m=0 

leading to 



m=0 



£ u m = f — l 1 1 r £ u m ) -l-Hn £ 



m=0 



m=Q 



(3) 



Now, in order to solve ([T| all one needs is to fix the components u m using Q. This 
can be done recursively, however, upon splitting first the action of Af on Y? n =o u m 
into sum of terms YZi=o^m in such a way that the first term Aq depends only on uq, 
which will be the first component of u to be fixed, and Ay depends only on (uq, u\), 
and A m depends only on (uo,uy, . . . ,u m ). With that said, the action of Af on u 
reads as 



A 7 £ u m — £ A m (uo,ui,. . .,u m ) . 



(4) 



m=0 



m=0 
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leading to 

OO / CO \ / CO \ 

£ u m =f-L- 1 ( £ 7? Mm ) -ZT 1 ( £ A m ) . (5) 

m=0 V m= / V m =0 / 

There is no unequivocal way to fix uq ||Wazwaz and El-Sayed |(|2001[)], however, 



a straightforward way to do it is to identify uq(x) with f(x). Similarly and in a 
straightforward way the other components are chosen to be 



U m + 1 



-L (Ru m )—L (A m ), m > 0, with uq = f. 



(6) 



If the action of N on u is represented by some function F(u), then A m (woi ^i ? ■ • • > ^m)> 
for m > 0, is expressed as a sum of m terms, each term includes a derivative of or- 
der k (1 < k < m) of F at mo, F^ k \uo), times a polynomial in (u\,U2, ■ ■ ■ ,u m ) and 
Ao = F(uo)- The first five so called Adomian "polynomials", Ao, . . . ,A4, available 
in the literature are 

Aq = F(uq), 
A\ = u\F'(uq) , 



A 2 = u 2 F'(u ) + ^yF"(i<o) , 



(7) 



Lf(3) 



A 3 = u 3 F'(u )+uiu 2 F"(u ) + j^F 

,2 



(wo) 



A 4 = m 4 F'(mo) + ( + "i«3 ]F"(»o) + 



~2\ 



^F^iu^ + P^i^U 



Adomian decomposition method (ADM) is known as a reliable mathematical tool 



for solving algebraic, differential or integral linear and nonlinear equations |E1 



Tawil, Bahnasawi, and Abdel-Naby (2004); 


Hosseini |(2006l; 


Kao, Jiang, and Yu 


(2005 ); Lai, Chen, and Hsu (2008 ); 


Wang, and Bajaj 


(2008); 


Wazwaz and El- 



Sayed ( |2001| ); Zhang (2005 1] to mention but a few. The method solves nonlinear 
equations appealing neither to linearization nor to perturbation approaches. ADM 
provides in many cases the exact solution, if it exists in closed form, or a semiana- 
lytical solution showing the main features of the exact one. However, a fastidious 
task in ADM is the evaluation of AP's which necessitates a big deal of calcula- 
tions. A first formal formula to evaluate AP's is that introduced by Adomian him- 
self |Adomian ( 1994)] followed by a series of some other alternative techniques 



developed in [Rach (1984); Seng, Abbaoui, and Cherruault (1996)] and some 
other papers aiming to reduce the size of calculations, however, they have fallen 
short of their aim. Based on the property that in each A m the sum of subscripts of 
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the components (wo, u\ u m ) is equal to m, a new algorithm has been developed 



to evaluate AP's without using Adomian formula [Wazwaz (2000)]. Besides the 
above mentioned property, AP's have some other interesting properties, never dis- 
cussed in the literature, which will be discussed and derived in the next section. 
The aim of this paper is to use these new properties as well as the above mentioned 
one to derive new reliable techniques for hand and computer evaluations of AP's. 
From the latter point of view a very simple and compact Mathematica program 
to compute any polynomial A m in its most general form as a function of F^\uq), 

1 < k < m, and of (u\,U2,.-.,u m ) is implemented. When the function F is given, 
the program derived below allows immediate evaluation of any individual poly- 
nomial A m without evaluating the polynomials with subscripts smaller than m and 
without using Adomian formula. Furthermore, the use of the new properties allows 
for elegant hand calculations of AP's. 

In section 2 we introduce the reduced polynomials Z m ^ and derive and discuss in 
detail their properties as well as the properties of A m . In section 3 we apply our 
results to derive general analytical expressions of A m for small values of m — k, 
which are valid for any F. In section 4 we discuss in detail the implementation 
of the Mathematica program and apply it to solve the problem of the pendulum in 
section 5. We conclude in section 6. 

2 Properties of Adomian polynomials 

The polynomials A m , m > 0, possess the following properties 

• Property 1. A m depends by construction only on the vector (u\,U2, ■ ■ ■ , u m ) 
and does not depend on u n with n> m; 

• Property 2. A m is the sum of m terms of the form 

m 

A m = ^ Zm,k(ui, U2,...,u m )F^\uo) (8) 

k=l 

where Z m ^ are called here the reduced polynomials and depend on (u\ , U2, ■ ■ ■ , 
u m ). We will see below that the properties of A m are encoded in Z m y, 

• Property 3. In each reduced polynomial Z mJ t, the components of the vec- 
tor (hi,M2, u m ) appear in such a way that the sum of their subscripts 
is equal to m [Wazwaz" ( 2000| )]. A consequence of this property is that 



each reduced polynomial Z m ^ depends only on few components of the vec- 
tor (u\,U2, ■ ■ ■ ,u m ). This is obvious since the sum of the subscripts of the 
components of the vector {u\ ,U2, ■ ■ ■ , u m ) exceeds m; 
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• Property 4. Another quit important property, curiously enough it has never 
been discussed in the literature to our knowledge, is that in each Z m ± the sum 
of the powers of the components u„, 1 < n < m, is equal to k (the order of 
the derivative in ([8])). Hence, both subscripts of Z m ^ have a direct meaning: 
one is the sum of the subscripts and the other is the sum of the powers of the 
components of the vector {u\,U2, ... ,u m ) which appear in the expression of 

The other properties are discussed below. Using the Taylor expansion for F(uq + e) 
in powers of e, one obtains 

(°° \ 00 F^Hu ) 

^) + [»iKW+£ (u l +u 2 + --) k . (9) 

A 

Now, substituting the multinomial formula 

rr n i 



(m+ii2 + "-) =kl 



ni +n 2 +---=k n/=l n j- 



into Eq. ([9]) and using the properties 3. and 4. of Z m ^ (m > 0), one obtains 



Zm,k — 7 


n i «2 


u m 


w n\ \ni\ ■ 


■n m \ 


with S 


= joi,. 


■ ■ i Km 



(10) 

In the above expression of Z mJ t, if some power nt = we have then uf = 1 and the 
corresponding component m, will not appear in the expression of Z m #. Now, since 
Z m> k depends only on few components of the vector (u\,U2, ■ ■ ■ ,u m ) as mentioned 



earlier, some of the powers in ( [10] ) must vanish. The summation in ( [10] ) being 
restricted to the set S, the powers are solutions to the system of diophantine 
equations 

«i+«2+«3H \-n m = k and n\ +2n 2 + 3«3 H \-mn m = m, (11) 



where we are interested to find solutions in N. Subtracting the first equation in ( [TT] ) 
from the second one, one obtains 

?i2 + 2«3 + 3«4 H h {m — k)n m -k+\ H h (wi — l)"m =m — k. (12) 

Since all the coefficients in (12) and the powers are positive integers or zero, 



Eq. ( [12] ), and consequently the system ( [TT] ), holds only if all the components n\, 
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i > m — k+ 1, vanish identically for all m and 1 < k < m. Said otherwise, all 
possible solutions to the system ( [TT] ) are such that: 



n m -k+2 — n m -k+3 — n m -k+4 



0, 



(13) 



meaning that Z m>k does never depend on (u m - k+2 ,u m - k+ 3, • • • ,u m ). We have thus 
established the proof that Z mk depends at most on (u\,u 2 , ■ ■ ■ ,u m -k+\), reducing 
Eq. (TlOb to 



l m-k+\ 



Sr 



n l \n 2 \---n m _ k+l \ 



with 



Sr 



(n u ...,n m - k+ i) £ N 



m-k+l 

"'~ k+1 \ £ 



m—k+l 
k and V in, 

1=1 



m 



(14) 



Said otherwise, the system ([TT]) is equivalent to the following reduced system of 
diophantine equations 

n\+n 2 -\ Vn m - k+ \=k and ni+2n 2 -\ h (m — k+ \)n m - k+ \ = m, (15) 



where Eq. ([T2j) becomes 

722 + 2«3 + 3«4 H h (m — 1 



m — k. 



(16) 



Does Z m £, as provided by ( [T4| ), depend on the whole vector (ui,u 2 , u m -k+i)^ 
We will shortly show that this is the case if 1 < k < m, however, if k = 1, Z m \ does 
not depend on the whole vector (u\ , u%, . . . , u m ); rather, it depends only on u m . 

Assume that 1 < k < m. It is straightforward to show that Z m _ k depends on u\ and 
linearly on u m - k+ \ for 1 < k < m. In fact, if we set n m - k+ \ = 1 in (fT6]> then the 



remaining powers in ([16]) must vanish and this leads using the first equation in ( [TT] ) 
to n\ = k — 1. Furthermore, the solution n\ = k — I and n m - k+ \ = 1 (all the other 
powers are zero) is the unique solution to ( j_5 1 for which n m - k+ \ / 0. Consequently, 
Z m k will always have the term 



u \ u m—k+l 



(17) 



in its expansion. For k = m, the term (17 1 reduces to «™/m!. To show that Z mJ t 
depends on u 2 and u m - k we first rewrite the numerator in ([T7]) as 



k—2 

U l U\U m - k+ \ 
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Now, rising the subscript of u\ by 1 and lowering the subscript of u m -k+\ by 1, in 
the above expression, while keeping u { unchanged, leads to 

k—2 k—2 

Notice that such a transformation does not affect both the sum of the powers and 
the sum of the subscripts. Hence, n\ = k — 2, n% = 1 and « m _* = 1 is a solution 
to ([B]) (or {□}) and 

k-2 

(.is; 



(*-2)l 



is a term in the expansion of Z m j c . Continuing so by rising one subscript and low- 



ering the other subscript in ( [T8] > while keeping u\ 2 unchanged, leads to 



k-2 



( ;_ 2 y ; ' 2) (2</</ max ) (i9) 

where / max = (m — k + 2)/2 if (m — /c) is even or / max = (m — k+ l)/2 if (m — k) 
is odd. Terms of the form ( fT9] > all appear in the expansion of Z m £. Now, as / runs 
from 2 to / max , the components m and u m _ k _^_ 2 ) span the whole vector (w2,«3, . . . , 
u m -k) meaning that Z m ^ depends on (w2,«3,..., « m -fc) an( l consequently, taking 
into account ( [TT] ), it depends explicitly on the whole vector (u\ ,U2,---, u m -k+i)- 

Assume k = 1 . The system of equations ( [T5] ) reduces to 

w i + n 2 H h «m = 1 and «i + 2«2 H h mn m = m . (20) 

Subtracting the second equation from the first one x m, one obtains 

(m — l)n\ + (m — 2)«2 H (wi — i)n; H \-1n m -%-\- n m -\ = 0, 

where the left-hand side is a sum of positive numbers with positive coefficients. 
Hence, the unique solution to the previous equation is n\ = n% = • • • = n m -\ = 0. 
Substituting into ( |20| ) leads to n m = 1. We have then shown that the unique solution 
to ([20} is n m = 1 and the other powers are zero, so Z MJj i depends only on u m and is 
given by 



-m.l 



(21) 



From the above discussion one can now write a more or less developed expression 
for A m (m > 0) 

A m = Z m> i(u m )F' (u ) H hZ m;fc (wi , « 2 , . . . , K TO _jfc+i (kq) 

H hZ mim _ 3 (Mi,M2,M3,«4)^ (m_3) («o) +Z m;m _ 2 («l,«2,W3)^ (m_2) ( M 0) 

+ Z m , m ^{u h u 2 )F^- l \u ) +Z nhm {u i )F (m \u ) . (22) 
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If one reads the formula ( |22| ) from the right to the left one sees that u\ appears for 
the first time in Z m ^ m , that U2 appears for the first time in Z m ,„,_i, that uj, appears for 
the first time in Z„,. m _2, u m -\ appears for the first time in Z m ^ and u m appears 
for the first time in Z m> \. 

One can still derive further properties of Z m % from ( fT6] > and ([15]) (or ( [T2] ) and ([TT])), 
however, for the analysis that will follow in the next section we will rely only on the 
properties discussed here. We can now apply the previous results to derive some 
expressions for Z m ^, m being fixed, which will serve as samples for the implemen- 
tation the Mathematica program in the next section. 



3 Applications and further properties 



According to (|22|) (or ( 14 1), Z mm depends only on u\ and to obtain the value of the 



power n\ we solve ( 15 ) or ( 1 1 1: n\ =m. Hence 



ml 



(23) 



Similarly, Z m m _i depends on u\ and ui. Solving ( [15] ) which reduces to: 
n\+ri2 = m— 1 and tii+2ti2 = 'w 



leads to «2 = 1 and n\ =m — 2. Hence 



m—2. 



Z mm —\ 



~U2 



(m-2)! 



(24) 



which is of the form < \TT} . Z m m _2 depends on (ui,U2,u^). The powers (n\, 112,113) 
satisfy the system of equations 



n\ +722 +«3 = m — 2 and n\ +2«2 + 3«3 = m . 



(25) 



which leads upon subtracting the first equation from the second one to 722 + 2«3 
2. This last equation has two possible solutions in N: («2 = 0,«3 = 1) or («2 



2,«3 = 0). We obtain the corresponding values of n\ from the first equation in ( |25] ): 
n 1 = m — 3 or n\ =m — 4, respectively. Hence 



„m-3„0„l 



,m-4„2,,0 



(m-3)!0!l! (m-4)!2!0! 



Mj W3 



,m-4„2 



+ 



(m-3)! (m-4)!2! 



(26) 
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Similarly, one obtains 



,m-4 



11 4 



,ni—5 



- , m,m—3 



Zm.m— 4 



Zm.m— 5 



(m-4)! 



+ 



M2M3 



*-6„3 



+ 



[m - 
,,m—6 



5)1 (m-6)!3!' 



U2U4 



.m-6 ,.2 



,m-1..2 



+ 



+ 



W 2 M 3 



+ ■ 



m-8,,4 



/iS 



(m-5)! (m-6)! (m-6)!2! (m-7)!2! (m-8)!4!' 



,m— 6. 



"6 



.hi— 7, 



+ 



W3W4 



.HI— 7, 



+ 



M2W5 



m-8,,2 



U2U4 



(m-6)!' (m-7)! ' (m - 7)! + (m - 8)!2! + (m - 8) !2! 



„m-9„3 



U2U3 



+ 



it 



m-10,.5 



- / m,\ 



(m-9)!3! (m-10)!5! ' 
u m , this has been derived in d21 



(27) 



We have also obtained the expression of Z m m _6 which is the sum of 1 1 terms. One 
can use the above formulas to evaluate the AP's A\ to Aj. For instance, if m = 7 
one obtains Z7 2 from the expression of Z„, m _5 dropping all negative powers of u\ 
since they are not in N. Hence, Z7 2 = «i«6 + M 3"4 + 1<2«5- Similarly, one obtains 
the other Z7 ^ (1 < k < 7) from the above formulas. 



A final conclusion, which can be read from ( |24| ), ( |26| ) and ( |27| ) (or from (|T4j)), is 
that each reduced polynomial Z m ^ is a sum of some number of terms, each term 
is a product of k factors and each factor represents one of the components of the 
vector («i,M2, • ■ • j Um-k+i)- For instance, Z7 2 depends on 7 — 2 + 1 = 6 functions: 
(u\,U2, U(,). Z-J2 is a sum of some terms (3 terms), each term is a product of 
k = 2 factors and each factor is one of the six functions (wi,t<2, ■ • • , «6) sucn that 
the sum of the powers is 2 and the sum of the subscripts is 7. This remark helps, in 
fact, improving hand calculations of A m . 



4 Implementation of the Mathematica program 

The Mathematica Reduce function helps solving diophantine equations. For in- 
stance, Eqs. (fT5]> for m = 4 and k = 2 write as 



(28) 



^rii = 2 and V in, = 4 . 



Mathematica solves (28 1 as follows: 



3 3 

Reduce[£n, == 2 && £in, == 4 && «i > && « 2 > && n 3 > 

i=l 1=1 

&& {ni,n2,«3} € Integers] 

(„i ==0&&?i2 ==2&& n 3 ==0)||(ni == 1 && « 2 ==0&& n 3 == 1) 
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There are two possible solutions and Z4.2 = (1*2/2!) + u\u 3 which is the coeffi- 
cient of F"(uq) in the fifth line of (|7]). In the above codes, the conditions n\ > 
&& «2 ^ && «3 > and {«i,«2j«3} G Integers depend on m and are not con- 
venient for large values of m. We replace them by the equivalent expressions 
3 3 

£Abs[n,] == £«,- and Table[«,-, {i, 3}] € Integers, respectively. These new ex- 

i=l i=l 

pressions are very important for the mathematica program we will develop in sub- 
section [43] where only m and F remain as free parameters. In the following two 



subsections, we provide two further applications of Reduce to solve < \15\ and ( TIT] ) 
for m = 4 and 1 < k < 4, respectively, and introduce and discuss in detail the codes 
needed for the Mathematica program. 

4.1 Solving Eqs. ([15]) for ra = 4 



In the following Mathematica codes we solve ( fT5] ) for m = 4 and 1 < k < 4, how- 
ever, the codes are valid for any value of m: changing m leads to new results. In the 
codes P represents the list of all possible power solutions. It is a list of four objects 
(elements) where each object is made of one or many sublists. For manipulating 
lists and sublists see pon 1 42000"] )] . 



m = 4; P = Table[0, {j,m}]; k = 0; 

m-k+l m-k+l 

While p = k + l)<m, P[[k}] = Reduce [ £ n, ==Jfc&& £ im==m 

i=l i=l 

m—k+l m—k+l 

&& Abs [«i] == £ "i && Table[«„ {i,m- k + l}} £ Integers]] 

i=l i=l 

P 

{n 4 == 1 &&n 3 == 0&&n 2 == 0&&7ii == 0, 

(m ==0&&«2 ==2&&n 3 ==0)||(m == 1&&« 2 ==0&&n 3 == 1), 
«1 == 2&&«2 == 1, 

m == 4} 

One sees P as a list of four objects, each object appears in one line. The objects 
P[[l]], P[[3]] and P[[4]] are each made of one sublist, while the object P[[2]] is 
made of two sublists separated by ||. The object P[[2]] corresponds to m = 4 and 
k = 2, which we know has two solutions. A drawback of solving ([15]) is that the 
different sublists P[[l]], P[[2, 1]], P[[2,2]], P[[3}} and P[[4}} have different lengths. 
For instance, Length[P[[l]]] = 4, Length[P[[2, 1]]] = 3 and Length[P[[3]]] = 2. We 
will see in the next subsection that solving ( fTT| ) will always lead to uniform sublists 
with the same length. 
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4.2 Solving Eqs. {FT} for m = 4 



In the following Mathematica codes we solve ( 1 1 1 for m = 4 and 1 < k < 4, how- 



ever, the codes are valid for any value of m: changing m leads to new results. 
„, = 4; p = Table[0, {j,m}] ;k = 0; 

in m 

While [(A: = k + \) <m, P[[k]] = Reduce == k && £/n, == m 

i=l 1=1 

m m 

&& ^Abs[«,-] == Y"i && Table[n„ {i,m}] € Integers]] 

i=l i=l 

{„ 4 == l &&n 3 == 0&&n 2 == 0&&7ii == 0, 
(„ 4 ==0&&« 3 ==0&&7i 2 ==2&&ni — — 0) 1 1 

(«4 == 0&&« 3 == 1 &&« 2 == 0&&«i == 1), 
M 4 == 0&&« 3 == 0&&722 == 1 &&«1 == 2, 
„ 4 == 0&&« 3 == 0&&712 == 0&&«l == 4} 

Notice that the sublists P[[l]], P[[2,l]], P[[2,2]], P[[3}] and P[[4]] have the same 
length 4 (= m). In order to be able to extract the values of Z^,* (1 < A: < 4) and 
evaluate A4 one needs to generate the corresponding list, U, of the functions u,. 
This is done by the code line U = P I. {n — *• «} which replaces 72 by w in the list of 
powers P. The list U has the same properties as the list P. 

Now, if the object U[[k]] of U (or the object P[[&]] of P) has only one sublist (these 
are the objects t/[[l]], U[[3]] and U[[4]] in our example) then 

m - k ~ njL, /»[[*, 7,2]]! ' 

and if the object t/[[&]] of U has many sublists (this is the object U[[2\] in our 
example), the number of sublists being Length [t/ [[&]]] leads to 

Length[f/[M]] n m =i TJ^iJ, l]]i>[[MJ,2]] 
Z m,k = 2- nm put. ; • on, • ( 3 °) 



-1 



TTJ =1 P[[M,./,2]]! 



4.3 TTie Mathematica program 

In order to distinguish between the case where U[[k]] has only one sublist and the 
case where it has many sublists we use the condition U[[k, 1,1,1]] ==0 which 
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holds if U[[k]] has only one sublist with m > 1. The Mathematica program is read- 
ily written for the case m = 7 which has been chosen for illustration: 

m = l;P = Table[0, {j,m}]; 
If[m> 1,4 = 0; 
While[(Jfc = Jfc+l)<m, 

mm mm 

P[[k]] = Reduce[ £ m ==£&&£ m, == m&& £ Abs[n f ] == £ m && 

i=i i=i i=i i=i 

Table[n/,{i,m}] € Integers]]; 
U = P/.{n^u}; 
A m ,i = 0;A m ,2 = 0; 
Do[If[Length[f/[[)fc, 1,1, 1]]] == 0, 

nU[[k,j,l)] p M 

A m ,i = A mA + ^ Derivative^] [F] [uq] , 

n(P[[*,i,2]]!) 

A m , 2 = A m , 2 + I —m Derivative^] [F] H], 

v w n«i,j,2]]o y 

{fe,l,m}]; 

A m = A m | + A ffli 2] 

If [m==l , A m = «i Derivative[l] [F] [wo]] 
If[m=0,A m = F[uq\] 



Aj = UtF'(uo) + (U3U4 + M2«5 +MlM6)F // (wo) 
,2„ „ „2 „2 



+ ^ + ^ +MiM2H4 + ^ f(3)(mo) 

— + — + — JF< >(«o) + + l4-J F< >o) 



^W+^'W (3D 



The constants that appear in the above expression are products of factorials: 5040 = 
7!, 12 = 3!2!, etc. In case of applications, the functions F and uq have to be defined 
earlier. 
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5 Application of the program: the pendulum 

A solution for the problem of the pendulum can be written down explicitly using 
elliptic functions, however, the pendulum is used here for mere illustration of the 
program since it is known to all workers. In the case where there are no friction and 
applied forces, the differential equation describing the motion of a pendulum is 

u tt + bsinu = 0, (32) 

where u(t) is the angular displacement and b is a positive "geometric" constant. 
The initial conditions are w(0) = a and u t (0) = 0, F(u) = sinu and the operator 
L _1 takes the form 

L- 1 [f(t)}= jUjy{p)}dp)dq. (33) 



Applying L 1 to both sides of ( 32 ) and using the prescribed initial conditions, one 
obtains 

u{t) = a-L~ l [bsm(u)}. (34) 
The problem of the pendulum is solved upon choosing 

uq = a and u m+ \ = —LT i [bA m ] , (m > 0) , (35) 



where A m are evaluated using the Mathematica program developed in subsection 4.3 
We have used ten polynomials A m : from Ao to A9 to evaluate eleven components of 
u(t): from uo(t) to uio(t). The results of calculations are shown below as produced 
by Mathematica 

10 

u = Collect[£ t/,,f] 
1=0 

bt 2 Sin\a] b 2 t 4 Cos\a]Sm\a] 



2 24 
6 /6 3 Sin[a] 6 3 Sin[3a]\ „/ b 4 Sin[2a] llb 4 Sin[4a} 



! ! 360 720 ) \ 5040 161280 

10 / b 5 Sin[a] 13^ 5 Sin[3a] 31^ 5 Sin[5a]\ 
1 V 45360 + 604800 3628800 J + 
12 /37ft 6 Sin[2a] _ 37ft 6 Sin[4a] 691Z7 6 Sin[6a] \ 
1 V 17107200 17107200 + 958003200 ) + 

/73Z7 7 Sin[a] b 7 Sin[3a} 9557^ 7 Sin[5a] 5461^ 7 Sin[7a] ^ 



340540200 3439800 43589145600 87178291200 
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16 f 313ft 8 Sin[2a] 37ft 8 Sin[4a] 9683ft 8 Sin[6a] 929569ft 8 Sin[8a] \ 
1 V 12573792000 + 1067489280 ~ 435891456000 + 167382319104000 J + 
18 / 47ft 9 Sin[a] t 381779ft 9 Sin[3a] 



20415732480 100037089152000 
816337ft 9 Sin[5a] 1441031ft 9 Sin[7a] 320229 lft 9 Sin [9a] \ ^ 



200074178304000 640237370572800 6402373705728000 y 
2Q f 1884343ft 1() Sin[2a] 540809ft 10 Sin[4a] 31786477ft 10 Sin[6a] 
1 V 63356823 12960000 ~ 1055947052160000 + 67580611338240000 ~ 
138706613ft 10 Sin[8a] 221930581ft 10 Sin[10a] \ 
608225502044160000 + 4865804016353280000 J ' 

The special case with a = n/2 reduces to 

71 ft? 2 b 3 t 6 b 5 t w lift 7 ? 14 211ft 9 ? 18 
U ^> ~ 2 ~ ~2 + 240 ~ 19200 + 13977600 ~ 16293888000 ' 

One may wonder whether the expression of u(t) for arbitrary a reduces to the linear 
case for small values of a. A series expansion of u(t,a) about the point a = and 
for an arbitrary value of t reads as produced by Mathematica 



Series[u,{a,0,3}] 

{ _bt^ bY_bW b 4 t % b 5 t lQ b 6 t 12 ft 7 ? 14 

~ ~2 + ~2A ~ 720 + 40320 ~ 3628800 + 479001600 ~ 87178291200 + 

ft 8 ? 16 ft 9 ? 18 b w t 20 \ 

20922789888000 ~ 6402373705728000 + 2432902008176640000 ) a + 

bt 2 b 2 t 4 5b 3 t 6 13ft 4 ? 8 1849ft 5 ? 10 4153ft 6 ? 12 149473ft 7 ? 14 
+ ■ 



12 36 864 15120 21772800 718502400 523069747200 
21019ft 8 ? 16 1100627ft 9 ? 18 5448101ft 10 ? 20 \ 3 

1961511552000 + 3492203839488000 ~ 729870602452992000 ) ° + ^ 



It is straightforward to check that the coefficient of a in the above expression 
is cos(\/ft?)- Hence, to the first order of approximation the solution is u(t) = 
acos(Vbt) which is the exact solution to the linear problem: u tt + bu = with 
the same initial conditions. 



6 Conclusion 

For a nonlinear operator N depending on one variable, it was shown that the reduced 
polynomial Z m ^ depends explicitly on m — k + 1 functions (mi,M2, . . • , u m -k+\) if 
k > 1 and that Z m .\ depends only on u m . The AP A m splits into a sum of m terms 
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where each term is the product of Z m ^, which is independent of N, with a derivative 
factor, FW(mo), depending on N. 

General expressions for Z m ^ have been derived for small values of m — k, hand 
calculations have been discussed briefly and a simple and compact Mathematica 
program has been implemented and applied. 

For the purpose of hand calculations and analytical studies, we complete the lists ([7]) 
and ( |3Tj ) providing general expressions for A m , which is not available in the litera- 
ture, by adding the polynomials A$, Ag, A% to A\$ using our program. 



A 5 

u 5 F'[u ] + (u 2 u 3 + u 1 ua)F"[u ] + ( ^ + ^ V (3) N + 4^ (4) M 



2 2 



A 6 



"l«2 , M l M 3\ F (4) N + ^ F (5) N + ; iiL F (6) W; 



24 J 720 



A* 

UgF'[u ] + f y + «3"5 + "2«6 + «1"7 ) f"[«o] 



+ 
+ 



( 2 2 2 \ 

+ — ^— + «1«3«4 + "1 «2«5 + I ™ ; [«0j 

4 2 2 2 2 3 \ 

9 UiU^Ut. U^U-j UiUtUa UtUs. \ „i A \ r , 



/ u\u\ u[u 2 U 3 u\ua \ „ f5U , 

+ \-n + ^r + ^4) F()[u ° ] 

fl44 I^UA (6) r i «f^2 (7) r i ^L F (8) [ 1 . 

V 48 120 / L J 720 L J 40320 L J 

A 9 

W)F' [uo] + («4«5 + «3«6 + «2«7 + "1 Mg)^" ["()] 

( 3 2 2 2 \ 

Mi M9W5 W1W7 \_,c<!i r i 

— + M 2 «3«4 + + — y + «1«3"5 + «1"2«6 + ~ y 1^ 1 [«0j 

( 3 2 2 2 2 3 \ 

MtW^ U]UlUq, U]U^Ua UiUt.Ua UiUlU*. UiUfr \ / A \ r 
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( u\u\ u\u]ut, u\u\ u\u 2 u A u\u 5 \ F ( 5 )r i 

V 24 4 12 6 24 ' ' 
f u\u\ u\u 2 m U \ U *\ F (6) h , 1 

V 36 24 120 y L 0J 

u\u\ u\llT,\ / 7 n m]m2 / 8 w , 



240 720 



2 



M1oF'[mo] + ^ y + M4«6 + U3UJ + M2«8 + U[U 9 j F"[ll ] 

+ I h — \-U 2 U 3 U 5 +U\U 4 U 5 + — h«lM3"6 + "l"2"7 + ^— )^ ; [«oJ 

( 2 2 3 3 22 
M1M3 «2 M 4 M 1 M 4 
— 1 ~ 1 ~ h M1M2W3M4 H — 
4 6 6 4 

+ Ml«2«5 + «i«3«5 + «f«2«6 + "|"7 \ p (4) ^ j 

2 2 2 6 / 

«2 UxU^U-i u\u 2 u\ u\u\u 4 u\ll3U4 U^U 2 Us u\uf,\ ,c\. 

T20 + ^ + ^ + ^ _ + ~6 _ + ~6 _ + ^a) F N 

, M^2 "i«2 M 3 «i«3 «f«2«4 M l "5 ^17(6)1- 1 

1 48 12 48 24 120 / [ "° J 

„4„3 „5„ „ „6„ \ / „6„2 „7 

1 144 



W3_ ^ \ ^\ 

120 720 y \ 1440 5040 / L 0J 



^ (9) N + ^dL^ (10) N. 



M i" 2 J7(9)r„j , M i 



40320 J 3628800 
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